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Abstract. It has recently been suggested that collapse of neutron stars induced 
by a phase transition to quark matter can be a considerable source of gravitational 
waves with kHz frequencies. We demonstrate that if about one percent of all neutron 
stars undergo this process, the resulting cosmological gravitational wave background 
would reach about 10~^° times the critical density. The background would peak 
at kHz frequencies and could have an observationally significant tail down to Hz 
frequencies. It would be comparable or higher than other astrophysical backgrounds, 
for example, from ordinary core collapse supernovae, from r-mode instabilities in 
rapidly rotating neutron stars, or from magnetars. The scenario is consistent with 
cosmological backgrounds in neutrinos and photons. 
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1. Introduction 

A phase transition from nuclear to quark matter within a neutron star (NS) with 
miUisecond scale rotation period can induce a collapse that releases up to ~ 10^^ erg in 
energy, which potentially relates this process to violent phenomena such as core collapse 
supernovae ^ or 7— ray bursts (GRBs) OlSlIll- Recent numerical simulations suggest 
that up to a few percent of this energy may be released as gravitational waves (GWs) 0. 

If a significant fraction of NSs undergo such a phase transition induced collapse, 
the energy released will contribute to cosmological backgrounds in photons, neutrinos, 
and GWs. In the present paper we estimate these backgrounds and point out that the 
stochastic GW background may be comparable to GW backgrounds from various other 
astrophysical sources. We specifically compare the power and statistical properties of 
this potential background with the ones from standard core collapse supernovae [Hj, NS- 
NS coalescence [7l|H], r-mode instabilities in NSs with miUisecond periods [Hllini, and 
from magnetars ^T], as well as from various processes in the early universe. We also 
find a significant constraint on the fraction of the total energy released in such phase 
transitions in form of MeV 7— rays. 

We use natural units with h = c = 1 throughout. 



2. The Source Mechanism 



We first recall that the energy radiated in GWs per frequency interval for an individual 
event at distance D, {dEg„/df){f), is related to the Fourier transform h{f) = 
I~oo e~^'^'^^^h{t) of the dimensionless strain amplitude h{t) by 

^(^'-TH^i^'"(^)^ (1) 

where Gn is Newton's constant. 

We model the strain in the phase transition scenario as 

hit) oc exp (-rt) sin f^t for t > , (2) 

which roughly reflects the form found in the simulations in Ref. 0. This leads to a 
Fourier transform 

Note that the total energy emitted in GWs from Eq. converges for this spectral 
form. We use /o = 2 kHz, F = 1/3 ms for the damping scale which is typical for the most 
rapidly rotating progenitors in Ref. |S|, and normalize such that ~ 2 x 10^^ erg, the 
more optimistic case from Ref. 

Most of the energy released in the phase transition is actually carried away in form 
of photons and neutrinos, whose total energy we will denote by E^ot ^ (1 ~5) x 10^^ erg. 
This is roughly the energy required to power a 7— ray burst. The gravitational wave 
energy is typically E^ ~ O.OSi^tot 
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Figure 1. Models for the star formation rate as functfon of redshift. The solid line is 
from Ref. jl5| which includes a putative PopIII component, whereas the dashed line 
respresents the Baldry & Glazebrook form fitted in Ref. 



3. Cosmological Event Rates 

The cosmic star-formation rate (SFR) and, as a consequence, the formation rate R{z) 
of NSs is reasonably well known at redshifts z < b [121 CH] • Especially for z < 1 it is 
strongly constrained by the Super-Kamiokande limit on the electron antineutrino flux 
from cosmological core collapse supernovae [131 lllj. In contrast, the SFR is poorly 
known for z > 5. We use the fits shown in Fig. 

We are interested in the number of events per unit mass. This number is given by 
an expression of the general form 

J dmm^im) ' 

where we neglect any redshift dependence. Here, S,{'m) is the mass function, usually 
taken to be of the Salpeter form ^(m) oc m~^-^^ between O.IMq and lOOM©, [Mmm, Mmax] 
is the progenitor mass range over which NSs form, and the integral in the denominator 
goes over all stellar masses. Further, x is the fraction of NSs eventually undergoing the 
phase transition under consideration. The event rate R{z) is then given by the product 
of the SFR and A. 

For NSs, Mmin ^ lOM©, M^^^ ~ 4OM0, so that A ~ 5 x lO^^x^©^- Note that 
A oc M~in^^ and could be significantly larger if NSs are also formed from progenitors 
with masses m < lOM©. The above value for A also roughly corresponds to today's 
rate of core collapse supernovae which is ^ 2 x 10^"^ Mpc^'^ yr^^ which is to be expected 
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given that most of the NSs are born in type II supernova events. Most of the hot 
young NSs are born with 10-20 ms rotation periods due to processes such as magnetic 
braking. We assume that x — 1% of the NSs undergo a phase transition. This is 
roughly the estimated fraction of NSs which are born with millisecond scale rotation 
periods [d HE] , required to trigger the phase transition. The same sub-population of 
NSs would be relevant for GW emission between ~ 100 Hz and ~ 2 kHz from r-mode 
instabilities in rapidly rotating NSs, in contrast to Ref. j^lCni which assumed x ~ 1- 

Another possible channel for phase transitions to occur in NSs is in low-mass X-ray 
binaries where the NS is spun up by accretion from the low mass companion |31 IT^ . 
The formation efficiency of such systems is hard to estimate, but typical numbers are 
X ~ 10~^ — 10~^ [20] 5 similar to the estimated fraction of millisecond NSs born in core 
collapse super novae. 

Even larger values for A could be motivated if these phase transitions are connected 
to short hard 7— ray bursts. A best estimate of ~ 10~^ Mpc~^ yr~^ for the comoving rate 
of faint short hard 7— ray bursts at zero redshift has recently been given Comparing 
with Fig. [T[ this corresponds to A ~ 5 x IO^'^Mq^, which would correspond to values 
for X comparable to unity. 

For the NS-NS coalescence and magnetar scenarios, A is actually given by the same 
expression Eq. (jH), with estimated fractions x = 0.3% |H| and x = 8% of all NSs 
contributing, respectively. 

We will neglect any time delay between star formation and NS formation because 
the lifetime ~ 10* yr of > 10 Mq NS progenitors is short compared to the Hubble time. 



4. The Gravitational Wave Background 



For simplicity, we assume that all NSs have identical emission characteristics. 

The energy density in GWs at frequency / per logarithmic frequency interval in 
units of the cosmic critical density pc = 3Hq/{8ttG-!<s) can be written as |221 



n,w(/) 



Pc Jo I + Z 



dt 



dz 



dE 

Jz {Jz 



(5) 



where R{z) is the event rate per comoving volume, fz = {l + z)f. The cosmological 
model enters with \dt/dz\ = [(1 -|- z)H{z)]~^ and, for a flat geometry, 

1/2 



Hz) 



Hn 



nM(i + zY + 



(6) 



We will use the parameters = 0.3, = 0.7, and Hq = ho 100 km s ^ Mpc ^ with 
ho = 0.72. 

For a stochastic GW background the dimensionless strain power fh{f) is related 
to the dimensionless GW energy density in Eq. (0) by 



ho ^lgw{f) 



fh{f) 



6.25 X 10 



35 



Hz 



fHf) 



(7) 
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Figure 2. Spectrum of GW background from Eq. 10) assuming about x — 1% of 
all NSs undergo a phase transition, releasing E^^ ^ 2 x IQ^"^ erg in GWs. The band 
represents the uncertainty related to the different SFRs shown in Fig. ^ The solid 
line is the maximal low-frequency tail due to anisotropic neutrino emission. For 
comparison, the dashed line is the maximal background from conventional type II 
supernovae discussed in Ref. |H] . Further, the horizontal dash-dotted lines represent a 
maximum version of the GW stochastic spectrum produced during slow-roll inflation 
assuming T / S = 0.3 for the ratio of the tensorial and scalar contributions to the 
cosmic microwave background radiation anisotropy and ±10"'^ for the running of the 
tensorial power-law index |^. The dotted lines show approximate sensitivities of 
the ground based interferometer LigoIII/EURO and LigoIII in correlation mode |24|. 
and of possible second generation space-based interferometers such as the Big Bang 
Observatory (BBC) (25| and the ultimate DECIGO as indicated. 



Fig. El shows the resulting GW background for the NS phase transition scenario, 
the band representing uncertainties due to the different SFRs from Fig. [TJ but not due 
to other quantities such as and x which were fixed to fiducial values, as indicated. 
Note that GW densities are proportional to A from Eq. and thus also to the fraction 
X of NSs subject to the phase transition. They are also proportional to the average total 
GW energy output -Egw As a result, the largest uncertainties come from the parameters 
X and i^gw, whereas uncertainties due to the SFR are moderate because the background 
is dominated by redshifts z <?> where the SFR is reasonably well known, see Fig. ^ 

The stochastic GW background from cosmological supernovae was studied recently 
in Ref. [HI, and its optimistic estimate is reproduced in Fig. |21 While the event rate 
R{z) used there is higher than in the present scenario (since supposedly only a percent 
fraction of core collapse supernovae likely give rise to NSs rotating rapidly enough to 
undergo phase transitions), the individual source signal is very different: In the phase 
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transition scenario, the signal from one event at kHz frequencies is much higher than 
the one due to convection in the simulations in Ref. [2Z1 and, as a consequence, the total 
energy emitted in GWs is much larger, ~ 10~^ Mq versus ~ 10~^ Mq. As a result, the 
background from phase transitions can be comparable or higher (above ~100 Hz) to the 
background from conventional type H supernovae. 

There could also be an enhanced low frequency tail if the strain h{t) converges to 
a non- vanishing constant for t ^ oo due to anisotropic neutrino emission ^Hl 1^ . For 
/ ^kHz, dEg^/df ^ 15Gti {E^q)"^ [0], where E^ < E^ot is the total energy emitted 
in neutrinos, and |g| < 1 is the average dimensionless quadrupole. Whereas E^ is 
about an order of magnitude smaller than i^j, ~ 3 x 10^^ erg in standard type H 
supernovae, the anisotropy could be much larger than the order percent anisotropy 
expected in hot NSs without phase transition j22j . This is because the NS would oscillate 
strongly and be highly deformed after the phase transition. In one of the simulations 
in Ref. 0, for example, the ratio of the polar to equatorial radius is ~ 0.7. As in core 
collapse supernovae, the neutrinos would likely be partially trapped and emitted from 
a neutrinosphere because the release of ~ 5 x 10^^ erg in internal energy would heat up 
the NS to ~ 10 - 20 MeV. 

We can parametrize the low-frequency tail as 

"-(^' - ' ' ' (57^) ^ (^TioW) (is) ■ 

The maximum of this tail, corresponding to the fiducial values in Eq. (jH)) is shown in 
Fig. 121 and reflects the uncertainty in the low-frequency signal: In the absence of a 
"GW memory effect", h{f) ^ const, see Eq. and thus VLg^{f) oc p for / ^ 0, see 
Eqs. and (0). In contrast, in the presence of burst with memory fh{f) ^ const and 
QgM) cx / for / ^ 0. 

Fig. El shows various other astrophysical backgrounds for comparison: The solid line 
shows a possible popcorn like signal from NS-NS coalescence jHI- We note that other 
work used a lower maximal frequency which would cut off the signal above ~ 200 Hz (Zj. 
The dotted line shows a stochastic background from r-mode instabilities emitting GWs 
between ~ 100 Hz and ~ 1.4 kHz. It has been corrected by the small fraction x ~ 1% of 
all NSs nowadays believed to be born with sufficiently small, millisecond scale rotation 
periods, in contrast to the original Refs. [011101 which assumed that a substantial fraction 
of all NSs are born with rotation periods close to maximal. Finally, the dashed line 
represents the magnetar scenario [llj. The NS phase transition scenario could thus 
produce a background which is comparable to or larger than these backgrounds both at 
kHz frequencies and below ~ 10 Hz. 

Finally, the background from NS phase transitions shown in Fig.Elwould potentially 
mask several possible GW relic backgrounds from the early universe, at least at 
frequencies between ~ 0.1 Hz and ~ 1 Hz, where it becomes close to gaussian and 
thus indistinguishable to other gaussian backgrounds, see next section. Such early 
universe backgrounds include the one from standard inflation whose maximum is shown 
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Figure 3. GW backgrounds from various other astrophysical mechanisms for 
comparison to Fig. ^ The solid hne represents that part of the NS-NS coalescence 
signal from Ref. which has a "popcorn" -character, similar to the phase transition 
scenario considered here, see below. The dotted line corresponds to GWs from r-mode 
instabilities emitted by the fastest rotating NSs 0^], assumed to constitute a fraction 
~ 1% of all NSs. The dashed line shows the magnetar scenario 



in Fig. 121 but also more speculative ones such as from quintessential inflation which 
predicts l]g„/ig ~ 6 x 10-i6(//Hz) above ~ IQ-^Hz jHO]. 

5. Statistical Properties of Gravitational Wave Background and 
Detectability 

The event rate as seen from Earth is 

/ dz^^—= dzR(z)- r^^, 9 

Jo 1 + z dz Jo ^ \l + z)H{z) ' ^ ' 

where dV/dz is the fractional volume element, the cosmic expansion rate at redshift z 
is given by Eq. ©, and r{z) is the comoving coordinate, dr = {1 + z)dt. 

The duty cycle can be estimated by multiplying the integrand of Eq. Q with 
(1 + z)tcoh[(l + -2)/], where tcoh(/) is the timescale over which the frequency / is emitted 
coherently. A probably more realistic estimate is obtained by weighting the resulting 
integral over redshift with the redshift dependent contribution to the GW signal given 
by the integrand in Eq. (jSj). 

Fig. 0] shows the duty cycle obtained by these two estimates for the NS phase 
transition scenario for which tcoh(/) — niax(l/r, 1//). If it is smaller than one, the 
duty cycle roughly corresponds to the variable ^ in Ref. The background is thus 
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Figure 4. The duty cycle as function of frequency for the NS phase transition scenario. 
The blue (lower) band includes weighting with the GW signal, whereas the dark (upper) 
band does not. Again, the bands represent the uncertainty related to the different SFRs 
shown in Fig. ^ 

of "popcorn" character. Note that the duty cycle is proportional to A. The minimal 
detectable GW energy density fig^ for a popcorn signal is comparable or, for duty 
factors smaller than a few percent, smaller by factors of a few compared to the minimal 
detectable background for gaussian noise Thus for A ~ 5 x 10~^Mq^, the signal 
would have a duty factor of a few percent and may thus be marginally detectable by 
third generation interferometers in correlation mode, see Fig. |21 This would require, 
however, that a fraction of order one of NSs would have to contribute, which is unlikely. 
Also note that microwave cavity detectors could reach sensitivities of order Qg^ ~ 10~^ 
at a few hundred Hz [321 ESI • 

The most advanced interferometers planned for the future such as GEO-HF j31] 
can reach sensitivities around f^^'^h{f) ~ 5 x 10-24 in the 1 - 10 kHz regime. 

One can easily deduce from Eqs. ((H) and (jHl) that for a source at distance D (small 
compared to the Hubble radius), the maximum of f^^'^h{f) occurs at / ~ 2kHz with 
p/^h{f) ~ 3 X 10-23 {D/W Mpc)-i Hz-^/2 ^^le phase transition scenario. Such events 
would thus be detectable out to ~ lOMpc, with a rate ~ 0.2 (x/0.01)(L'/10Mpc)2 yr-^ 

At / ~ kHz the signal from core collapse supernovae has a duty cycle slightly higher 
than the one of the phase transition signal if x ~ 1% of the NSs contribute to the phase 
transition signal. This is because the supernova rate is then roughly a factor 100 higher, 
whereas the coherence time tcoh(/) — 1// ^^e incoherent process of convection is 
slightly shorter than the one of the oscillations right after the phase transition, for which 
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For both the r-mode instabihty OCni and magnetar [TT] scenarios, tcoh(/) ^ 1//, 
and the signal is Gaussian at all observable frequencies. The NS-NS coalescence signal 
consists of a gaussian, a popcorn, and a shot noise component jH]. A fraction of 
X ~ 0.003 of all NSs are in suitable binaries to contribute to the coalescence signal. 
As a consequence, the background has statistical properties similar to the one for NS 
phase transitions. The total power from NS-NS coalescence shown in Fig. 01 can be 
comparable or somewhat higher than for the NS phase transition scenario because the 
gravitational energy released is higher, Eg„ ~ 5 x 10^^ erg jH]. 

6. Electromagnetic and Neutrino Emissions 

The energy going into GWs during the phase transitions considered in Ref. [5] represents 
only a fraction of a few percent of the total energy released £'tot- Fig. |21 shows that an 
energy Eg^ ~ 2 x 10^^ erg released per event in GWs and a number of events per unit 
mass X 5 X 10"^Mq^ corresponds to figw ~ 10^^°, quite independent of the poorly 
known high redshift evolution of the SFR. By simple scaling, this implies an energy 
density 



is released in form of photons and/or neutrinos. Photons and neutrinos are expected 
to be emitted on timescales > 10 s ^ l/F than GWs, a factor > 3 x 10^ larger than 
GWs which are emitted on a timescale ~ l/F ~ 3 ms. Fig. EJthen implies that these 
backgrounds are quasi-continuous for A > 5 x 10"^ Mq^, provided the emissions are not 
strongly beamed. 

This emission can certainly not be predominantly in form of MeV 7— rays because 
the cosmological diffuse background of photons above an MeV has an energy density 
^Mev^o — 4 X 10^^°. Put another way, the fraction /mbV that can be released in these 
events is constrained by 



Since the energy density in lower energy backgrounds are larger, no significant 
constraint on the fraction going into photons of energy below ~ 100 keV result. 

Obviously, if only a very small fraction of the SFR contributes to these events, i.e. 
if A is small, and/or the total energy released i^tot is small, this constraint becomes weak 
or disappears. However, the fiducial values for A and i^tot used here are not unrealistic, 
and imply significant constraints on emission into MeV 7— rays. 

If phase transitions are to be associated with GRBs, a fraction 




(10) 




(11) 




(12) 
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released as ~ MeV 7— rays would in fact be sufficient because GRB emissions correspond 
to a (non-continuous) photon background with energy density ^grb ~ 2 x 10~^^. Since 
~ 10^ GRBs are visible per year, this would imply a beaming factor 

~ 3 X 10=^ f ^] . (13) 

The corresponding Lorentz factor would be roughly the square root of this. 

7. Conclusions 

We have estimated the cosmological background of gravitational waves for the scenario 
where a fraction x ~ 1% of all neutron stars are born as millisecond pulsars and 
undergo a phase transition to a quark star, releasing ~ 10^^ erg in gravitational waves, as 
suggested by recent simulations. In terms of the corresponding event rate per unit stellar 
mass. A, this background would have a duty factor of ~ 5 x 10~^ A/(5 x 10~^Mq^ 



above ~ 50 Hz and ~ 0.03 



A/(5 X 10-5 M-i; 



(//Hz) ^ below ~ 50 Hz. However, most 



of the energy flux is at kHz frequencies, a factor 10-100 above the maximum sensitivity 
of ground based interferometers. Using statistics speciahzed to detecting popcorn type 
noise with collocated aligned detectors at the technological limit, the background may 
be marginally detectable around 100 Hz for A ~ 5 x 10^'^ ^0 ^ which is only possible if 
the majority of all neutron stars would be born with millisecond periods and undergo a 
phase transition. In this unrealistic case the background would have duty factors of a 
few percent. 

The total energy in gravitational waves can be, however, quite substantial even for 
moderate fractions x of neutron stars contributing. In units of the critical density one 
obtains flgwh^ ~ 10"^° [£'gw/(2 x 10^^ erg)] (x/0.01). In addition, since neutrinos tend 
to be emitted very anisotropically in the form of jets, there could be a low-frequency tail 
of order fig^(/) < 10^^^ A/ (5 x 10^^ ^^o^) (//^z), likely larger than for ordinary type 
II supernovae. A low-frequency tail of that size would mask the maximal gravitational 
background in ordinary inflation models between ~ 0.1 Hz and ~ 1 Hz where it would 
also be gaussian. At these frequencies this tail could be detectable by future space based 
interferometer projects such as BBO or DECIGO. 

The gravitational wave background from phase transitions in rapidly rotating 
newly born neutron stars can thus be comparable or higher than other astrophysical 
backgrounds. Detection of such a signal could provide insight into the nature of compact 
objects. 

If most of the energy released in the phase transition is in form of neutrinos, 
the resulting diffuse MeV neutrino flux would constitute ~ 10~^(x/0.01) of the flux 
due to standard type II core collapse supernovae, and would thus be consistent with 
existing upper limits ^3]. The fraction of the energy released in form of MeV 7— rays 
is constrained by observed diffuse backgrounds to be less than ~ 10%(x/0.01)"^. A 
fraction ~ 10-^(x/0.01)-^ released into the band of 100 keV — 1 MeV 7— rays would 
suffice to account for 7— ray bursts. 
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